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Abstract 

Nuclear pairing is studied both in atomic nuclei and in neutron-star crusts in the unified frame- 
work of the energy-density functional theory using generalized Skyrme functionals complemented 
with a local pairing functional obtained from many-body calculations in homogeneous nuclear 
matter using realistic forces. 
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I. INTRODUCTION 



Thepossibility of pairing in atomic nuclei was first studied by Bohr, Mottelson and 
Pines jl] and by Belyaev {^J only one year after the publication of the theory of supercon- 
ductivity by Bardeen, Cooper and Schrieffer (BCS) p]. Meanwhile, Bogoliubov developed a 
microscopic theory of superfluidity and superconductivity and explored its consequences for 
nuclear matter [4|. In 1959, Migdal speculated that the interior of neutron stars might be su- 
perfluid J5J and this scenario was further investigated by Ginzburg and Kirzhnits in 1964 jg]. 
Soon after the discovery of the first pulsars, the observation of frequency glitches followed 
by very long relaxation times of the order of months provided strong evidence of nuclear 
superfluidity [7]. Pulsar glitches are believed to be related to the dynamics of the neutron 
superfluid permeating the inner layers of the solid neutron star crust 8|. Superfluidity plays 
also a predominant role in neutron-star cooling (see Page in this volume). 

The pairing phenomenon in both finite systems like atomic nuclei and in infinite nu- 
clear matter can be consistently described in the framework of the energy- density functional 
(EDF) theory (see Dobaczewski and Nazarewicz in this volume). This theory, which has 
been historically formulated in terms of effective interactions in the context of self-consistent 
mean-field methods, has been very successful in describing the structure and the dynam- 
ics of a wide range of nuclei js]. These interactions have been also commonly applied to 
the modeling of neutron-star interiors. Actually no sooner did Skyrme [lo| introduce his 
eponymous effective interaction than Cameron llj applied it to calculate the structure of 
neutron stars. By showing that their maximum mass was significantly higher than the Chan- 
drasekhar mass limit, his work brought support to the scenario of neutron star formation 
from the catastrophic gravitational core-collapse of massive stars during type II supernova 
explosions, as proposed much earlier by Baade and Zwicky 12]. 



II. NUCLEAR ENERGY DENSITY FUNCTIONAL THEORY IN A NUTSHELL 

Assuming time-reversal symmetry, the ground-state energy E is supposed to depend on 
(i) the nucleon density (denoting the spin states by a = ±1 and q = n or p for neutron or 
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proton, respectively), 

Pl( r ) = Pq( r ^^^) > (!) 

cr=±l 

(ii) the kinetic-energy density (in units of h 2 /2M q where M q is the nucleon mass), 

)= E / dV5(r-r')V-V'p g (r,a;r',cr), (2) 



T q (r) 



cr=±l 

(iii) the spin-current density, 

J q (r) = -i I dV 6 ( r - r ')Vp g (r, a; r', a 1 ) x cv CT (3) 

<t,ct'=±1 ^ 

and (iv) the abnormal density, 

Pi( r ) = Pq(r,v,r,a) , (4) 

cr=±l 

where oVo-' denotes the Pauli spin matrices. In turn the normal and abnormal density 
matrices, p(r, cr; r', a') and p(r,<7;r' ,<j') respectively, can be expressed as 

p 9 (r,cr;r>') = (r,*)^ (r', </)* (5) 



and 



p,(r,a;rV) = -X;^if (r.^SPC/yr, (6) 

<(ff) 



where (r, cr) and 4i'( r i (T ) are the two components of the quasiparticle (q.p.) wave- 



function. Here, as throughout this paper pure nucleon states are being assumed^ 



;he more 



general formalism involving neutron-proton mixing has been developed in Ref. 13] . 

Minimizing the total energy E with respect to ip[f(r, cr) and tp^(r,a) under the con- 
straints of fixed particle numbers leads to the Hartree-Fock-Bogoliubov (HFB) equations 1 
(see Dobaczewski and Nazarewicz in this volume) 

/ h q (r) aa , A q (r)5 a A U^(r,a') 
V U,(r)cw -h q (rU] \4t(r,v') 



E i + X q \ faiiir,*)' 

^-A g J USW) 



(7) 



1 These equations are also called Bogoliubov-de Gennes equations in condensed matter physics. 



where X q are Lagrange multipliers. The single-particle (s.p.) Hamiltonian h q (r) aa > is given 
by 

h q {r) c , a = -V ■ B q (r)V 5 aa > + U q {r)5 aa , - iW q {r) • V x a a , a (8) 
with the s.p. fields defined by the functional derivatives 

Using a local pairing EDF of the form 

WO = ~ E ^ 9 [Pn(r),p p (r)]p 9 (r) 2 . (10) 
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q=n,p 



the pairing field is given by 



rp -r 

A,M = ^ = ^"[ft.M.ft.MlftW . (li) 



fl r 



Expressions for these fields can be found for instance in Refs. |l5j. Eqs. fl7|)-f lTT]) are still 
valid at finite temperatures, but Eqs. ([5]) and ((6]) will have to be replaced by 

(12) 

and 

p>, a;r', a') = £(2/, - 1)4? (r, a)^\r', a'f , (13) 

i(s) 

where f\ are the q.p. occupation probabilities given by (setting the Boltzmann constant 
k B = l) 

fl= 1 + exp^/T) ■ (M) 

III. SKYRME FUNCTIONALS 

The nuclear EDF that we consider here is of the Skyrme type 9J, i.e., 

E = E kin + -Ecoui + E Sky + E pair , (15) 



where E^ n is the kinetic energy of the normalization volume, -Ecoui is the Coulomb energy 
(dropping the exchange part in order to simulate neglected effects such as Coulomb correla- 
tions, c harg e-symmetry breaking of the nuclear forces and vacuum polarization as discussed 
in Ref. [la]), -Esky is the Skyrme energy and -E pa ; r is the nuclear pairing energy. 

Historically the Skyrme energy was obtained from the Hartree-Fock approximation using 



an effective interaction of the form 



9] 



= * (1 + xoPMra) + + [p% S( rij ) + Sinjtfj] 

+t 2 {l + x 2 P a )^ 1 Pi j .5(r ij )pi j + ^t 3 (l + x 3 P a )p{r) a <J(r y ) 

+#o(ffi + &j) • Pij x 5(rij) Pij , (16) 
n 

where = r, — r,-, r = (fj + r 3 -)/2, = —ih(Vi — Vj)/2 is the relative momentum, P CT is 
the two-body spin-exchange operator. Likewise, the pairing energy can be obtained from a 
Skyrme-like effective interaction given by 

<f = ^(1 " P CT )^ 9 [Pn(r),p p (r)] 5(r y ) . (17) 

Because of the zero-range of the pairing force, a cutoff has to be used in the gap equations 
in order to avoid divergences (for a review of the various prescriptions, see for instance 
Ref. (17J). 

The density dependence of the pairing strength v" Kq [p n ,p p ] remains very poorly known. 
It has been usually assumed that it depends only on the isoscalar density p — p n + p p and 
has often been parametrized as [3] 

^v,p P ] = ^(i-%(£) 9 ) , (is) 

where p is the nuclear saturation density while VA, r\ q and a q are adjustable parameters. 
The superscript A on V^ q serves as a reminder that the pairing strength depends very strongly 
on the cutoff. In principle changing the cutoff modifies also the other parameters but the 
effects are generally found to be small. Effective interactions with rj q = (rj q = 1) have 
been traditionally refered as volume (surface) pairing. The standard prescription to fix the 



parameters is to adjust the value of the pairing strength V^ q to the average gap in 120 Sn 19]. 
However this does not allow an unambiguous determination of the remaining parameters r\ q 
and a q . Systematic studies of nuclei seem to favor a so-called mixed pairing with r] q ~ 0.5 



and l/2<a q <l |2o|,|2!| 



The parameters of the Skyrme EDF are usually determined so as to reproduce a set of 
nuclear data selected according to a specific purpose. The non-uniqueness of the fitting 
procedure has led to a large number of different parametrizations. Some of them may yield 
very different predictions when applied outside the domain where they were fitted. This 
situation is particularly unsatisfactory for nuclear astrophysical applications which require 
the knowledge of nuclear masses for nuclei so neutron rich that there is no hope of measuring 
them in the foreseeable future; such nuclei play a vital role in the r-process of nucleosyn- 
thesis 22] and are also found in the outer crust of neutron stars [23|. Extrapolations far 



beyond the neutron drip line are required for the description of the inner crust of neutron 



stars 24| where nuclear clusters are embedded in a sea of unbound neutrons, which are ex- 
pected to become superfluid at low enough temperatures. Even though the crust represents 
only ~ 1% of the neutron-star mass, it is intimately related to many observed astrophysi- 
cal phenomena j^j]. The need for more reliable extrapolations of these nuclear EDFs has 
motivated recent efforts to construct non-empirical effective interactions and more generally 
microscopic nuclear EDFs [26]. Unfortunately such ab initio nuclear EDFs have not yet 
been developed to the point where they can reproduce nuclear data with the same degree of 
accuracy as do phenomenological EDFs, which can now fit essentially all the nuclear mass 
data with rms deviations lower than 0.6 MeV [27]. 

IV. PAIRING IN INFINITE HOMOGENEOUS NUCLEAR MATTER 

In infinite homogeneous matter the HFB equations ([7]) can be readily solved. The q.p. 
energies are given by 



E ( k q) = v (4 9) - \y + ^ , (19) 

in terms of the s.p. energies 

ef = B q k 2 + U q . (20) 

The q.p. wavefunctions reduce to 

4i(r, a) = Mr, a) , V$(r, a) = fc (r, a) , (21) 

with 

1 / M ^ \ 1/2 

(22) 




-(<?) 



1/2 



V2 



E 



(<?) 



(23) 



and 0fc(r, a) is given by 



1 



exp (ifc • r) x((r) 



(24) 



where x( cr ) is the Pauli spinor and V is the normalization volume. The uniform pairing 
obeys the well-known isotropic BCS gap equations (see for instance Appendix B of Ref. 

E(e) 



he 



A q = - l -v^[ PniPp ] A q ^de^tanh. 



2T 



(25) 



where g(e) is the density of s.p. states (per unit energy) given by 

*) = 55F- (26) 

Dq 

The subscript A is to indicate that the integral has to be regularized by introducing a cutoff. 
We include here all s.p. states whose energy lies below X q + E\ where E\ is a pairing energy 
cutoff. 



In the weak-coupling approximation, i.e. A q <C X q and A q <C e\, it was shown in Ref. 28 1 
that the pairing gap at T = is approximately given by 



A q (0) = 2X q exp 



g{X q )v*i 



exp 



2 V A„ 



where 



A(x) = ln(16x) + 2Vl + x - 2 In (l + Vl + x^j - 4 . 
The Lagrange multiplier X q is approximately given by the Fermi energy 



Xq — B q k Fq , 



(27) 



(28) 



(29) 



with k-pq = (37r 2 p g ) 1 / 3 . Note that these expressions were obtained by going beyond the 
usual "weak-coupling approximation" in which the density of s.p. states is taken as a 
constant. Even though this provides a good approximation in the case of conventional BCS 
superconductivity [3j, it is less accurate in the nuclear context because many more states 
are involved in the pairing mechanism. The temperature-dependence of the pairing gap can 



be very well represented by 



29) 



A„(T<T c )^A g (0)Wl 



(30) 



with 5 ~ 3.23 and the critical temperature is given by 

T, 



A,(0)^Uo.57A,(0). 



7T 



(31) 



Phenomenological pairing functional whose parameters have been fitted to nuclei gen- 
erally yield unrealistic pairing gaps in homogeneous nuclear matter [14 



30]. Given the 



uncertainties regarding pairing correlations in nuclei, Garrido et al. {31] proposed to deter- 
mine the parameters of the pairing strength in Eq. (fT8l) by fitting the 1 5'o pairing gaps in 
infinite symmetric nuclear matter (SNM) as obtained by the realistic Paris potential in the 
BCS approximation. The pairing interaction between two nucleons inside a nucleus is thus 
assumed to be locally the same as the pairing interaction between two nucleons in infinite 
uniform matter. Even though the coupling to surface vibrations is expected to contribute 
to pairing 32J (see also Kamerdzhiev and Avdeenkov in this volume), a local pairing theory 



seems a reasonable first step (see Bulgac in this volume). The main difficulty of this ap 



proach is that because of the highly non-linear character of pairing correlations |33l . |34j it is 
very difficult to guess an appropriate functional form for v 7Tq [p n ,p p \ (see also Lombardo in 
this volume). For instance, Margueron et al. [35j have recently shown that the parametric 
form (ITS]) has to be generalized in order to reproduce the 1 Sq pairing gaps in both SNM and 
pure neutron matter (NeuM) as obtained from Brueckner calculations 



Alternatively, Eq. f l27|) can be inverted to obtain the analytic expression of the pairing 
strength in terms of a given pairing gap function A q (p n , p p ) 



V* q [pn,Pp] 



8tt 2 B. 



2r3/2 



Iq(Pn,Pp) 



with 



A 



£a 
A, 



1/ j 



(32) 



(33) 



The value of the pai ring cutoff e\ is not completely arbitrary and could be fixed as follows. 
It has been argued 3JJ, |37| that in the limit p — > 0, the pairing strength should coincide 
with the bare force in the 1 So channel, which in turn is determined by the experimental 1 So 
nucleon-nucleon phase shifts. However at very low densities, the pairing strength is simply 
given by 

r 2 / *2 \3/2 



v nq [p^ 0] 



4vr 2 



h 2 



el \2M Q 



(34) 
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FIG. 1: Root mean square deviation a between theoretical and experimental masses for some 260 
known spherical or quasi-spherical nuclei for different values of the pairing cut-off energies. See 
text for more details. 

The optimal value of the cutoff is thus found to be E\ ~ 7 — 8 MeV (note that E\ is half 
the cutoff used in Ref. (37|). On the other hand, such a low cutoff may not be optimal in 
applications to finite nuclei. As seen in Fig. [TJ the root mean square (rms) deviation obtained 
with respect to about 260 known masses of (quasi-) spherical nuclei is found to oscillate as 
a function of the cutoff energy with clear minima lying around 7, 16 and 24 MeV. In this 
example, the initial EDF has the same characteristics as BSk21 with a pairing functional 



constrained on nuclear matter properties 



271 ]. as described in the next section. Note that for 



each value of the cut-off, the Skyrme interaction parameters were re-adjusted to minimize the 
rms deviation. In particular, we found systematically that global fits to nuclear masses favor 
£ a ~ 16 MeV, a value which we adopted. Similar results were obtained when considering a 
traditional ^-pairing force with or without the Bulgac-Yu regularization [381 ]. 

V. PAIRING IN NUCLEI 



We have recently constructed a family of three nuclear EDFs, BSkl9, BSk20 and 
BSk21 27|] based on Skyrme forces that are generalized in the sense that they contain density- 



dependent generalizations of the usual t\ and t<i terms, respectively [15| . The neutron-pairing 
functional was obtained from Eqs. (!32|) - (l33|) using the 1 S pairing gaps both in SNM mat- 
ter and NeuM, as obtained from Brueckner calculations including medium polarization ef- 
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FIG. 2: Neutron pairing strength as a function of the nucleon density p = p n + p p and isospin 
asymmetry 77 = (p n - p p )/p. 



fects [36]. The resulting pairing strength is shown in Fig. [2j The proton-pairing functional 
had the same form but we allowed its strength to be different from the neutron-pairing 
strength in order to take account of Coulomb effects not included in the above nuclear mat- 
ter calculations. Because of our neglect of polarization effects in odd nuclei due to our use 
of the equal-filling approximation |39|, we also allowed each of these strengths to depend 
on whether there is an even or odd number of nucleons of the charge type in question 2 . 
These extra degrees of freedom were taken into account by multiplying the neutron-pairing 
functional v nq [p n , p p ], as determined by the nuclear-matter calculations that we have just 
described, with renormalizing factors f^, where /+, f~ and /~ are free, density-independent 
parameters to be included in the mass fit, with /+ set equal to 1. 

The remaining parameters of the functionals were determined primarily by fitting essen- 
tially all measured nuclear masses. For this it was necessary to add two phenomenological 
corrections to the HFB ground-state energy: (i) a Wigner energy (which contributes signifi- 
cantly only for light nuclei or nuclei with N close to Z) and (ii) a correction for the spurious 
rotational and vibrational collective energy. However, in fitting the mass data we simul- 
taneously constrained the functional to fit the zero-temperature equation of state (EOS) 
of NeuM, as determined by three different many-body calculations using realistic two- and 
three-nucleon forces. Finally, we imposed on these EDFs the supplementary constraints 

2 Note that the odd nucleon will nevertheless contribute to the time-even fields. 
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of i) eliminating all spurious spin-isospin instabilities in nuclear matter both at zero and 



finite temperatures, at all densities found in neutron stars and supernova cores [151. 1271. |40|. 
ii) obtaining a qualitatively realistic distribution of the potential energy among the four 
spin-isospin channels in nuclear matter, iii) ensuring that the isovector effective mass is 
smaller than the isoscalar effective mass, as indicated by both experiments and many-body 
calculations and iv) restricting the incompressibility coefficient of SNM to lie in the range 
K v = 240 ± 10 MeV, as inferred from breathing-mode measurements. 

The introduction of the unconventional Skyrme terms allowed us to satisfy all these 
constraints and at the same time fit the 2149 measured masses of nuclei with N and Z > 



8 given in the 2003 Atomic Mass Evaluation (AME) [41| with an rms deviation as low 
as 0.58 MeV for all three EDFs. Incidentally, our EDFs are found to be consistent with 
measurements of the high-density pressure of SNM deduced from heavy-ion collisions 42], 



even though they were not directly fitted to the EOS of SNM. 

VI. PAIRING IN NEUTRON-STAR CRUSTS 

Because of the precision fit to masses and the constraints on both the EOS and the 

n 

1 Sq pairing gaps in NeuM, our recently developed EDFs BSkl9, BSk20 and BSk21 [27| are 
particularly well-suited for describing the inner crust of neutron stars. 

The HFB equations flj[J) have been already solved in neutron-star crusts using the so- 
called Wigner-Seitz (W-S) approximation according to which the crust is divided into a 



set of independent spheres centered around each lattice site [43M45| (see also Barranco et 
al. in this volume). However this way of implementing the HFB method can only be 
reliably applied in the shallowest region of the inner crust where nuclear clusters are very 



far apart 461 ] . In order to investigate pairing correlations in the densest part of the crust, 
we have applied the band theory of solids, which takes into account both short- and long- 
range correlations 0|- The band theory relies on the assumption that the solid crust can 
be treated as a perfect crystal, which is a reasonable approximation for cold non-accreting 



neutron stars 



251 ] . According to the Floquet-Bloch theorem, the q.p. wave function must 



obey the following boundary conditions 48j for any lattice translation vectors I 

4ak(r +i,<r)= exp(ifc • t) ^(r, a) (35) 
11 



where a is the band index (principal quantum number) and k the Bloch wave vector. This 
formalism naturally incorporates the local rotational symmetries around the nuclear clusters 
as well as the translational symmetry of the lattice, thus describing consistently both clusters 
and superfluid neutrons. Note that this formalism also includes infinite homogeneous matter 
as the limiting case of an "empty" lattice. The band theory therefore allows for a unified 
treatment of all regions of a neutron star. 

In the deep layers of the inner crust of a neutron star, where spatial inhomogeneities 
are small, further simplifications can be made. In the decoupling approximation, the q.p. 
wavefunction is expressed in terms of the s.p. wavefunctions ip^l as 

1>£k<r, a) = U<$ <p$(r, a) , <^ = Vff ^ (r, a) . (36) 
The HFB equations can then be readily solved, and the q.p. energies are given by 



^ = V( £ S-A ? ) 2 +(A2) 2 (37) 
where e^t are the s.p. energies and are solutions of the anisotropic multi-band BCS 



gap equations 



47] 



\'"> _^ V f d3fc/ V® f^tanh^ (3S) 



with 



/ dV^nPn(r),p P (r)]|^(r)| 2 |^,(r)| 2 . (39) 
./ws 

The subscript WS indicates that the integral has to be taken inside the W-S cell. Finally 
the amplitudes of the q.p. wavefunction are given by 




1/2 

(40) 




VS^1-^J • (41) 
We have solved Eq. (1381) for neutrons in the deep region of neutron-star crusts as described 



in Ref. 



471 ] using our latest BSk21 EDF (which is strongly favored by the most recent 



atomic mass data while being also consistent with what is now known about neutron-star 
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FIG. 3: Neutron superfluidity in neutron-star crust at average baryon density p = 0.06 fm~ 3 with 
BSk21. The critical temperature is found to be equal to T c = 1.1 MeV. Left panel: pairing gaps 
vs single-particle energies at T = 0. The arrow indicates the position of the Fermi energy. Right 
panel: pairing field in the W-S cell obtained from the HFB equations (solid line) and from the 
LDA (dashed line) for two different temperatures. 



masses 



49||). Due to the presence of the nuclear clusters, neutrons belonging to different 
bands and having different Bloch wave vectors feel different pairing interactions thus leading 
to a dispersion of the neutron pairing gaps of a few hundred keV around the Fermi 
level, as shown in Fig. [3j The critical temperature is found to be very weakly dependent on 
the cutoff, as can be seen in Table fl] (note that e\ was varied while using the same Skyrme 
functional BSk21). 

Because the neutron superffuid coherence length is much larger than the size of the 
clusters, proximity effects are very strong. As a result, pairing correlations are substantially 
enhanced inside clusters while they are reduced in the intersticial region, leading to a smooth 
spatial variation of the pairing field. The local density approximation (LDA), whereby the 
neutron pairing field A n (r) is assumed to be locally the same as that in uniform nuclear 
matter for the neutron density p n (r) and proton density p p {r), strongly overestimates the 
spatial variations of the pairing field. The discrepancies are particularly large inside clusters 
where the LDA incorrectly predicts a quenching of pairing correlations, especially for tem- 
peratures close to the critical temperature as illustrated in Fig. [31 This analysis shows that 
a consistent treatment of both unbound neutrons and nucleons bound in clusters is essential 
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TABLE I: Cut-off dependence of the critical temperature of neutron superfluidity in neutron star 
crusts at density p = 0.06 fm -3 with BSk21. See text for further details. 



£A [MeV] 


2 


4 


8 


16 


32 


T c [MeV] 


1.10 


1.06 


1.07 


1.11 


1.16 



for a realistic description of pairing correlations in neutron-star crusts. 

Despite the absence of viscous drag at T = 0, the solid crust can still resist the flow 
of the neutron superfluid due to non-dissipative entrainment effects. These effects have 
been systematically studied in all regions of the inner crust of a cold non-accreting neutron 
star [50|. In particular, it has been found that in some layers of the inner crust, almost all 
neutrons are entrained by clusters. These results suggest that a revision of the interpretation 
of many observable astrophysical phenomena like pulsar glitches might be necessary 51]. 



VII. CONCLUSIONS 



The nuclear EDF theory opens the way to a unified description of the nuclear pairing 
phenomenon in various systems, from atomic nuclei to neutron stars. 

The Brussels-Montreal EDFs based on generalized Skyrme EDFs supplemented by a 
microscopic local pairing EDF yield an excellent fit to essentially all experimental nuclear 
mass data with rms deviations falling below 0.6 MeV, while reproducing at the same time 
many-body calculations in infinite homogeneous nuclear matter using realistic forces. For 
this reason, these EDFs are particularly well-suited for studying pairing correlations in the 
inner crust of neutron stars, where nuclear clusters are expected to coexist with a neutron 
superfluid. 

Despite these successes, a number of open issues like for instance neutron-proton pairing 
or the contribution of surface vibrations to pairing call for more elaborate pairing EDFs. 
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